Folate pathway polymorphisms and neuropsychological impairment after leukemia therapy

Neuropsychological impairment occurs in 20%-40% of childhood acute lymphoblastic leukemia (ALL) survivors, possibly mediated by folate depletion following methotrexate chemotherapy. We evaluated the relationship between two folate pathway polymorphisms and neuropsychological impairment after childhood ALL chemotherapy. Eighty-six childhood ALL survivors were recruited between 2004-2007 at Texas Children's Hospital after exclusion for central nervous system leukemia, cranial irradiation, and age<1 year at diagnosis. Neuropsychological evaluation at a median of 5.3 years off therapy included a parental questionnaire and the following child performance measures: Trail Making Tests A and B, Grooved Pegboard Test Dominant-Hand and Nondominant-Hand, and Digit Span subtest. We performed genotyping for polymorphisms in two folate pathway genes: reduced folate carrier (RFC1 80G>A, rs1051266) and dihydrofolate reductase (DHFR Intron-1 19bp deletion). Fisher exact test, logistic regression, Student's t-test, and ANOVA were used to compare neuropsychological test scores by genotype, using a dominant model to group genotypes. In univariate analysis, survivors with cumulative methotrexate exposure ≥9000 mg/m2 had an increased risk of attention disorder (OR=6.2, 95% CI 1.2 – 31.3), compared to survivors with methotrexate exposure <9000 mg/m2. On average, female survivors scored 8.5 points higher than males on the Digit Span subtest, a test of working memory (p=0.02). The RFC1 80G>A and DHFR Intron-1 deletion polymorphisms were not related to attention disorder or impairment on tests of attention, processing speed, fine motor speed, or memory. These data imply a strong relationship between methotrexate dose intensity and impairment in attention after childhood ALL therapy. We did not find an association between the RFC1 80G>A or DHFR Intron-1 deletion polymorphisms and long-term neuropsychological impairment in childhood ALL survivors.^

GSTM1 and GSTT1 null polymorphisms and risk of salivary gland carcinoma.

Glutathione S-transferase (GST) genes detoxify and metabolize carcinogens, including oxygen free radicals which may contribute to salivary gland carcinogenesis. This cancer center-based case-control association study included 166 patients with incident salivary gland carcinoma (SGC) and 511 cancer-free controls. We performed multiplex polymerase chain reaction-based polymorphism genotyping assays for GSTM1 and GSTT1 null genotypes. Odds ratios (ORs) and 95% confidence intervals (CIs) were calculated with multivariable logistic regression analyses adjusted for age, sex, ethnicity, tobacco use, family history of cancer, alcohol use and radiation exposure. In our results, 27.7% of the SGC cases and 20.6% of the controls were null for the GSTT1 (P = 0.054), and 53.0% of the SGC cases and 50.9% of the controls were null for the GSTM1 (P = 0.633). The results of the adjusted multivariale regression analysis suggested that having GSTT1 null genotype was associated with a significantly increased risk for SGC (odds ratio 1.5, 95% confidence interval 1.0-2.3). Additionally, 13.9% of the SGC cases but only 8.4% of the controls were null for both genes and the results of the adjusted multivariable regression analysis suggested that having both null genotypes was significantly associated with an approximately 2-fold increased risk for SGC (odds ratio 1.9, 95% confidence interval 1.0-3.5). The presence of GSTT1 null genotype and the simultaneous presence of GSTM1 and GSTT1 null genotypes appear associated with significantly increased SGC risk. These findings warrant further study with larger sample sizes.

Prognostic significance of xCT polymorphisms and expression in patients with advanced pancreatic cancer treated with chemotherapy

The plasma membrane xc- cystine/glutamate transporter mediates cellular uptake of cystine in exchange for intracellular glutamate and is highly expressed by pancreatic cancer cells. The xCT gene, encoding the cystine-specific xCT protein subunit of xc-, is important in regulating intracellular glutathione (GSH) levels, critical for cancer cell protection against oxidative stress, tumor growth and resistance to chemotherapeutic agents including platinum. We examined 4 single nucleotide polymorphisms (SNPs) of the xCT gene in 269 advanced pancreatic cancer patients who received first line gemcitabine with or without cisplatin or oxaliplatin. Genotyping was performed using Taqman real-time PCR assays. A statistically significant correlation was noted between the 3' untranslated region (UTR) xCT SNP rs7674870 and overall survival (OS): Median survival time (MST) was 10.9 and 13.6 months, respectively, for the TT and TC/CC genotypes (p = 0.027). Stratified analysis showed the genotype effect was significant in patients receiving gemcitabine in combination with platinum therapy (n = 145): MST was 10.5 versus 14.1 months for the TT and TC/CC genotypes, respectively (p = 0.013). The 3' UTR xCT SNP rs7674870 may correlate with OS in pancreatic cancer patients receiving gemcitabine and platinum combination therapy. Paraffin-embedded core and surgical biopsy tumor specimens from 98 patients with metastatic pancreatic adenocarcinoma were analyzed by immunohistochemistry using an xCT specific antibody. xCT protein IHC expression scores were analyzed in relation to overall survival in 86 patients and genotype in 12 patients and no statistically significant association was found between the level of xCT IHC expression score and overall survival (p = 0.514). When xCT expression was analyzed in terms of treatment response, no statistically significant associations could be determined (p = 0.908). These data suggest that polymorphic variants of xCT may have predictive value, and that the xc- transporter may represent an important target for therapy in pancreatic cancer.

Characterization of glucocorticoid receptor (NR3C1) gene polymorphisms and a family-based association study with hypertension

Hypertension is usually defined as having values of systolic blood pressure ≥140 mmHg, diastolic blood pressure ≥90 mmHg. Hypertension is one of the main adverse effects of glucocorticoid on the cardiovascular system. Glucocorticoids are essential hormones, secreted from adrenal glands in circadian fashion. Glucocorticoid's effect on blood pressure is conveyed by the glucocorticoid receptor (NR3C1), an omnipresent nuclear transcription factor. Although polymorphisms in this gene have long been implicated to be a causal factor for cardiovascular diseases such as hypertension, no study has yet thoroughly interrogated the gene's polymorphisms for their effect on blood pressure levels. Therefore, I have first resequenced ∼30 kb of the gene, encompassing all exons, promoter regions, 5'/3' UTRs as well as at least 1.5 kb of the gene's flanking regions from 114 chromosome 5 monosomic cell lines, comprised of three major American ethnic groups—European American, African American and Mexican American. I observed 115 polymorphisms and 14 common molecularly phased haplotypes. A subset of markers was chosen for genotyping study populations of GENOA (Genetic Epidemiology Network of Atherosclerosis; 1022 non-Hispanic whites, 1228 African Americans and 954 Mexican Americans). Since these study populations include sibships, the family-based association test was performed on 4 blood pressure-related quantitative variables—pulse, systolic blood pressure, diastolic blood pressure and mean arterial pressure. Using these analyses, multiple correlated SNPs are significantly protective against high systolic blood pressure in non-Hispanic whites, which includes rsb198, a SNP formerly associated with beneficial body compositions. Haplotype association analysis also supports this finding and all p-values remained significant after permutation tests. I therefore conclude that multiple correlated SNPs on the gene may confer protection against high blood pressure in non-Hispanic whites. ^

Role of obesity, weight gain, physical activity, and polymorphisms in the mTOR pathway and the risk of renal cell carcinoma

The interplay between obesity, physical activity, weight gain and genetic variants in mTOR pathway have not been studied in renal cell carcinoma (RCC). We examined the associations between obesity, weight gain, physical activity and RCC risk. We also analyzed whether genetic variants in the mTOR pathway could modify the association. Incident renal cell carcinoma cases and healthy controls were recruited from the University of Texas MD Anderson Cancer Center in Houston, Texas. Cases and controls were frequency-matched by age (±5 years), ethnicity, sex, and county of residence. Epidemiologic data were collected via in-person interview. A total of 577 cases and 593 healthy controls (all white) were included. One hundred ninety-two (192) SNPs from 22 genes were available and their genotyping data were extracted from previous genome-wide association studies. Logistic regression and regression spline were performed to obtain odds ratios. Obesity at age 20, 40, and 3 years prior to diagnosis/recruitment, and moderate and large weight gain from age 20 to 40 were each significantly associated with increased RCC risk. Low physical activity was associated with a 4.08-fold (95% CI: 2.92-5.70) increased risk. Five single nucleotide polymorphisms (SNPs) were significantly associated with RCC risk and their cumulative effect increased the risk by up to 72% (95% CI: 1.20-2.46). Strata specific effects for weight change and genotyping cumulative groups were observed. However, no interaction was suggested by our study. In conclusion, energy balance related risk factors and genetic variants in the mTOR pathway may jointly influence susceptibility to RCC. ^

Polymorphisms, linkage and mapping of four enzyme loci in the fish genus Xiphophorus (Poeciliidae).

Electrophoretic variants at four additional enzyme loci--two esterases (Est-2, Est-3), retinal lactate dehydrogenase (LDH-1) and mannose phosphate isomerase (MPI)--among three species and four subspecies of fish of the genus Xiphophorus were observed. Electrophoretic patterns in F1 hybrid heterozygotes confirmed the monomeric structures of MPI and the esterase and the tetrametric structure of LDH in these fishes. Variant alleles of all four loci displayed normal Mendelian segregation in backcross and F2 hybrids. Recombination data from backcross hybrids mapped with Haldane's mapping function indicate the four loci to be linked as Est-2--0.43--Est3--0.26--LDH-1--0.19--MPI. Significant interference was detected and apparently concentrated in the Est-3 to MPI region. No significant sex-specific differences in recombination were observed. This group (designated linkage group II) was shown to assort independently from the three loci of linkage group I (adenosine deaminase, glucose-6-phosphate dehydrogenase, and 6-phosphogluconate dehydrogenase) and from glyceraldehyde-3-phosphate dehydrogenase and two isocitrate dehydrogenase loci. Evidence for conservation of the linkage group, at least in part, in other vertebrate species is presented.

Polymorphisms of the DNA repair gene MGMT and risk and progression of head and neck cancer.

Methylating agents are involved in carcinogenesis, and the DNA repair protein O(6)-methylguanine-DNA methyltransferase (MGMT) removes methyl group from O(6)-methylguanine. Genetic variation in DNA repair genes has been shown to contribute to susceptibility to squamous cell carcinoma of the head and neck (SCCHN). We hypothesize that MGMT polymorphisms are associated with risk of SCCHN. In a hospital-based case-control study of 721 patients with SCCHN and 1234 cancer-free controls frequency-matched by age, sex and ethnicity, we genotyped four MGMT polymorphisms, two in exon 3, 16195C>T and 16286C>T and two in the promoter region, 45996G>T and 46346C>A. We found that none of these polymorphisms alone had a significant effect on risk of SCCHN. However, when these four polymorphisms were evaluated together by the number of putative risk genotypes (i.e. 16195CC, 16286CC, 45996GT+TT, and 46346CA+AA), a statistically significantly increased risk of SCCHN was associated with the combined genotypes with three to four risk genotypes, compared with those with zero to two risk genotypes (adjusted odds ratio (OR)=1.27; 95% confidence interval (CI)=1.05-1.53). This increased risk was also more pronounced among young subjects (OR=1.81; 95% CI=1.11-2.96), men (OR=1.24; 95% CI=1.00-1.55), ever smokers (OR=1.25; 95%=1.01-1.56), ever drinkers (OR=1.29; 95% CI=1.04-1.60), patients with oropharyngeal cancer (OR=1.45; 95% CI=1.12-1.87), and oropharyngeal cancer with regional lymph node metastasis (OR=1.52; 95% CI=1.16-1.89). In conclusion, our results suggest that any one of MGMT variants may not have a substantial effect on SCCHN risk, but a joint effect of several MGMT variants may contribute to risk and progression of SCCHN, particularly for oropharyngeal cancer, in non-Hispanic whites.

Association of retinoic acid receptor genes with meningomyelocele.

BACKGROUND: Neural tube defects (NTDs) occur in as many as 0.5-2 per 1000 live births in the United States. One of the most common and severe neural tube defects is meningomyelocele (MM) resulting from failed closure of the caudal end of the neural tube. MM has been induced by retinoic acid teratogenicity in rodent models. We hypothesized that genetic variants influencing retinoic acid (RA) induction via retinoic acid receptors (RARs) may be associated with risk for MM. METHODS: We analyzed 47 single nucleotide polymorphisms (SNPs) that span across the three retinoic acid receptor genes using the SNPlex genotyping platform. Our cohort consisted of 610 MM families. RESULTS: One variant in the RARA gene (rs12051734), three variants in the RARB gene (rs6799734, rs12630816, rs17016462), and a single variant in the RARG gene (rs3741434) were found to be statistically significant at p < 0.05. CONCLUSION: RAR genes were associated with risk for MM. For all associated SNPs, the rare allele conferred a protective effect for MM susceptibility.

Association of folate receptor (FOLR1, FOLR2, FOLR3) and reduced folate carrier (SLC19A1) genes with meningomyelocele.

BACKGROUND: Meningomyelocele (MM) results from lack of closure of the neural tube during embryologic development. Periconceptional folic acid supplementation is a modifier of MM risk in humans, leading toan interest in the folate transport genes as potential candidates for association to MM. METHODS: This study used the SNPlex Genotyping (ABI, Foster City, CA) platform to genotype 20 single polymorphic variants across the folate receptor genes (FOLR1, FOLR2, FOLR3) and the folate carrier gene (SLC19A1) to assess their association to MM. The study population included 329 trio and 281 duo families. Only cases with MM were included. Genetic association was assessed using the transmission disequilibrium test in PLINK. RESULTS: A variant in the FOLR2 gene (rs13908), three linked variants in the FOLR3 gene (rs7925545, rs7926875, rs7926987), and two variants in the SLC19A1 gene (rs1888530 and rs3788200) were statistically significant for association to MM in our population. CONCLUSION: This study involved the analyses of selected single nucleotide polymorphisms across the folate receptor genes and the folate carrier gene in a large population sample. It provided evidence that the rare alleles of specific single nucleotide polymorphisms within these genes appear to be statistically significant for association to MM in the patient population that was tested.

Influence of host interleukin-10 polymorphisms on development of traveler's diarrhea due to heat-labile enterotoxin-producing Escherichia coli in travelers from the United States who are visiting Mexico.

Up to 60% of U.S. visitors to Mexico develop traveler's diarrhea (TD), mostly due to enterotoxigenic Escherichia coli (ETEC) strains that produce heat-labile (LT) and/or heat-stable (ST) enterotoxins. Distinct single-nucleotide polymorphisms (SNPs) within the interleukin-10 (IL-10) promoter have been associated with high, intermediate, or low production of IL-10. We conducted a prospective study to investigate the association of SNPs in the IL-10 promoter and the occurrence of TD in ETEC LT-exposed travelers. Sera from U.S. travelers to Mexico collected on arrival and departure were studied for ETEC LT seroconversion by using cholera toxin as the antigen. Pyrosequencing was performed to genotype IL-10 SNPs. Stools from subjects who developed diarrhea were also studied for other enteropathogens. One hundred twenty-one of 569 (21.3%) travelers seroconverted to ETEC LT, and among them 75 (62%) developed diarrhea. Symptomatic seroconversion was more commonly seen in subjects who carried a genotype producing high levels of IL-10; it was seen in 83% of subjects with the GG genotype versus 54% of subjects with the AA genotype at IL-10 gene position -1082 (P, 0.02), in 71% of those with the CC genotype versus 33% of those with the TT genotype at position -819 (P, 0.005), and in 71% of those with the CC genotype versus 38% of those with the AA genotype at position -592 (P, 0.02). Travelers with the GCC haplotype were more likely to have symptomatic seroconversion than those with the ATA haplotype (71% versus 38%; P, 0.002). Travelers genetically predisposed to produce high levels of IL-10 were more likely to experience symptomatic ETEC TD.

Theoretical and experimental studies of linkage disequilibrium in human populations

Linkage disequilibrium (LD) is defined as the nonrandom association of alleles at two or more loci in a population and may be a useful tool in a diverse array of applications including disease gene mapping, elucidating the demographic history of populations, and testing hypotheses of human evolution. However, the successful application of LD-based approaches to pertinent genetic questions is hampered by a lack of understanding about the forces that mediate the genome-wide distribution of LD within and between human populations. Delineating the genomic patterns of LD is a complex task that will require interdisciplinary research that transcends traditional scientific boundaries. The research presented in this dissertation is predicated upon the need for interdisciplinary studies and both theoretical and experimental projects were pursued. In the theoretical studies, I have investigated the effect of genotyping errors and SNP identification strategies on estimates of LD. The primary importance of these two chapters is that they provide important insights and guidance for the design of future empirical LD studies. Furthermore, I analyzed the allele frequency distribution of 26,530 single nucleotide polymorphisms (SNPs) in three populations and generated the first-generation natural selection map of the human genome, which will be an important resource for explaining and understanding genomic patterns of LD. Finally, in the experimental study, I describe a novel and simple, low-cost, and high-throughput SNP genotyping method. The theoretical analyses and experimental tools developed in this dissertation will facilitate a more complete understanding of patterns of LD in human populations. ^

Nonlinear tests for genetic studies of complex disease

Linkage and association studies are major analytical tools to search for susceptibility genes for complex diseases. With the availability of large collection of single nucleotide polymorphisms (SNPs) and the rapid progresses for high throughput genotyping technologies, together with the ambitious goals of the International HapMap Project, genetic markers covering the whole genome will be available for genome-wide linkage and association studies. In order not to inflate the type I error rate in performing genome-wide linkage and association studies, multiple adjustment for the significant level for each independent linkage and/or association test is required, and this has led to the suggestion of genome-wide significant cut-off as low as 5 × 10 −7. Almost no linkage and/or association study can meet such a stringent threshold by the standard statistical methods. Developing new statistics with high power is urgently needed to tackle this problem. This dissertation proposes and explores a class of novel test statistics that can be used in both population-based and family-based genetic data by employing a completely new strategy, which uses nonlinear transformation of the sample means to construct test statistics for linkage and association studies. Extensive simulation studies are used to illustrate the properties of the nonlinear test statistics. Power calculations are performed using both analytical and empirical methods. Finally, real data sets are analyzed with the nonlinear test statistics. Results show that the nonlinear test statistics have correct type I error rates, and most of the studied nonlinear test statistics have higher power than the standard chi-square test. This dissertation introduces a new idea to design novel test statistics with high power and might open new ways to mapping susceptibility genes for complex diseases. ^

Genome-wide association study of LDL-cholesterol over ApoB ratio

Genome-Wide Association Study analytical (GWAS) methods were applied in a large biracial sample of individuals to investigate variation across the genome for its association with a surrogate low-density lipoprotein (LDL) particle size phenotype, the ratio of LDL-cholesterol level over ApoB level. Genotyping was performed on the Affymetrix 6.0 GeneChip with approximately one million single nucleotide polymorphisms (SNPs). The ratio of LDL cholesterol to ApoB was calculated, and association tests used multivariable linear regression analysis with an additive genetic model after adjustment for the covariates sex, age and BMI. Association tests were performed separately in African Americans and Caucasians. There were 9,562 qualified individuals in the Caucasian group and 3,015 qualified individuals in the African American group. Overall, in Caucasians two statistically significant loci were identified as being associated with the ratio of LDL-cholesterol over ApoB: rs10488699 (p<5 x10-8, 11q23.3 near BUD13) and the SNP rs964184 (p<5 x10-8 11q23.3 near ZNF259). We also found rs12286037 ((p<4x10-7) (11q23.3) near APOA5/A4/C3/A1 with suggestive associate in the Caucasian sample. In exploratory analyses, a difference in the pattern of association between individuals taking and not taking LDL-cholesterol lowering medications was observed. Individuals who were not taking medications had smaller p-value than those taking medication. In the African-American group, there were no significant (p<5x10-8) or suggestive associations (p<4x10-7) with the ratio of LDL-cholesterol over ApoB after adjusting for age, BMI, and sex and comparing individuals with and without LDL-cholesterol lowering medication. Conclusions: There were significant and suggestive associations between SNP genotype and the ratio of LDL-cholesterol to ApoB in Caucasians, but these associations may be modified by medication treatment.^

Genetic Predictors of Clinical Outcomes in Non-Small Cell Lung Cancer Patients

Lung cancer is the leading cause of cancer-related mortality in the US. Emerging evidence has shown that host genetic factors can interact with environmental exposures to influence patient susceptibility to the diseases as well as clinical outcomes, such as survival and recurrence. We aimed to identify genetic prognostic markers for non-small cell lung cancer (NSCLC), a major (85%) subtype of lung cancer, and also in other subgroups. With the fast evolution of genotyping technology, genetic association studies have went through candidate gene approach, to pathway-based approach, to the genome wide association study (GWAS). Even in the era of GWAS, pathway-based approach has its own advantages on studying cancer clinical outcomes: it is cost-effective, requiring a smaller sample size than GWAS easier to identify a validation population and explore gene-gene interactions. In the current study, we adopted pathway-based approach focusing on two critical pathways - miRNA and inflammation pathways. MicroRNAs (miRNA) post-transcriptionally regulate around 30% of human genes. Polymorphisms within miRNA processing pathways and binding sites may influence patients’ prognosis through altered gene regulation. Inflammation plays an important role in cancer initiation and progression, and also has shown to impact patients’ clinical outcomes. We first evaluated 240 single nucleotide polymorphisms (SNPs) in miRNA biogenesis genes and predicted binding sites in NSCLC patients to determine associations with clinical outcomes in early-stage (stage I and II) and late-stage (stage III and IV) lung cancer patients, respectively. First, in 535 early-stage patients, after correcting multiple comparisons, FZD4:rs713065 (hazard ratio [HR]:0.46, 95% confidence interval [CI]:0.32-0.65) showed a significant inverse association with survival in early stage surgery-only patients. SP1:rs17695156 (HR:2.22, 95% CI:1.44-3.41) and DROSHA:rs6886834 (HR:6.38, 95% CI:2.49-16.31) conferred increased risk of progression in the all patients and surgery-only populations, respectively. FAS:rs2234978 was significantly associated with improved survival in all patients (HR:0.59, 95% CI:0.44-0.77) and in the surgery plus chemotherapy populations (HR:0.19, 95% CI:0.07-0.46).. Functional genomics analysis demonstrated that this variant creates a miR-651 binding site resulting in altered miRNA regulation of FAS, providing biological plausibility for the observed association. We then analyzed these associations in 598 late-stage patients. After multiple comparison corrections, no SNPs remained significant in the late stage group, while the top SNP NAT1:rs15561 (HR=1.98, 96%CI=1.32-2.94) conferred a significantly increased risk of death in the chemotherapy subgroup. To test the hypothesis that genetic variants in the inflammation-related pathways may be associated with survival in NSCLC patients, we first conducted a three-stage study. In the discovery phase, we investigated a comprehensive panel of 11,930 inflammation-related SNPs in three independent lung cancer populations. A missense SNP (rs2071554) in HLA-DOB was significantly associated with poor survival in the discovery population (HR: 1.46, 95% CI: 1.02-2.09), internal validation population (HR: 1.51, 95% CI: 1.02-2.25), and external validation (HR: 1.52, 95% CI: 1.01-2.29) population. Rs2900420 in KLRK1 was significantly associated with a reduced risk for death in the discovery (HR: 0.76, 95% CI: 0.60-0.96) and internal validation (HR: 0.77, 95% CI: 0.61-0.99) populations, and the association reached borderline significance in the external validation population (HR: 0.80, 95% CI: 0.63-1.02). We also evaluated these inflammation-related SNPs in NSCLC patients in never smokers. Lung cancer in never smokers has been increasingly recognized as distinct disease from that in ever-smokers. A two-stage study was performed using a discovery population from MD Anderson (411 patients) and a validation population from Mayo Clinic (311 patients). Three SNPs (IL17RA:rs879576, BMP8A:rs698141, and STK:rs290229) that were significantly associated with survival were validated (pCD74:rs1056400 and CD38:rs10805347) were borderline significant (p=0.08) in the Mayo Clinic population. In the combined analysis, IL17RA:rs879576 resulted in a 40% reduction in the risk for death (p=4.1 × 10-5 [p=0.61, heterogeneity test]). We also validated a survival tree created in MD Anderson population in the Mayo Clinic population. In conclusion, our results provided strong evidence that genetic variations in specific pathways that examined (miRNA and inflammation pathways) influenced clinical outcomes in NSCLC patients, and with further functional studies, the novel loci have potential to be translated into clinical use.

Association of p21/p27 polymorphisms with risk of HPV16-associated oral squamous cell carcinoma

The increasing incidence of oral squamous cell carcinoma (OSCC) among young adults has been associated with sexually transmitted infection of human papillomavirus (HPV), particularly HPV16. Given the roles of p21 (WAF1/Cip1/CDKN1A) and p27 (Kip1/CDKNIB) in cell-cycle regulation and of HPV16 E6 and E7 oncoproteins in p53 degradation and pRb inactivation, the effect of HPV16 L1 seropositivity and three putatively functional single-nucleotide polymorphisms (SNPs) of p21 (p21 C70T and p21 C98A) and p27 (p27 T109G), individually and in combination, on the risk of OSCC was evaluated in a hospital-based case-control study of 327 cases and 401 cancer-free controls who were frequency-matched on age, gender and smoking status. Individuals with HPV16 L1 seropositivity had an overall 3-fold increased risk of having OSCC than those with HPV16 seronegativity. The increased risk of HPV16-associated OSCC was particularly found among younger people (aged ≤ 50 years), males, never smokers, never drinkers and oropharynx cancer patients. None of three p21 and p27 polymorphisms alone was significantly associated with risk of OSCC. Individuals with variant genotypes for both p21 polymorphisms were more likely to have OSCC than individuals with wild-type genotypes or variant genotypes for either one of the p21 polymorphisms (adjusted OR, 1.4; 95% CI, 0.9-2.1). There was a borderline significant or significant interaction between the p21 C70T, combined p21 and combined p21/p27 genotypes and HPV16 L1 seropositivity on risk of OSCC. The three studied p21 and p27 polymorphisms, individually or in combination, did not appear to have an effect on HPV16-related clinical outcomes (overall and disease-free survival and tumor recurrence). Despite the fact that the exact biological mechanism remains to be explored, these findings suggest possible involvement of p21variants, particularly the p21 C70T variant genotypes (CT/TT), in the etiology of HPV16-associated OPSCC. Further large and functional studies are required to validate the findings.^